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O ABSTRACT 

We present and interpret new X-ray data for M33SNR21, the brightest X-ray 

supernova remnant (SNR) in M33. The SNR is in seen projection against (and 

appears to be interacting with) the bright Hll region NGC592. Data for this 

■ source were obtained as part of the Chandra ACIS Survey of M33 (ChASeM33) 
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Very Large Project. The nearly on-axis Chandra data resolve the SNR into a 
~ 5" diameter (20 pc at our assumed M33 distance of 817 ± 58 kpc) slightly 
elliptical shell. The shell is brighter in the east, which suggests that it is en- 
countering higher density material in that direction. The optical emission is 
coextensive with the X-ray shell in the north, but extends well beyond the X-ray 
rim in the southwest. Modeling with X-ray spectrum with an absorbed sedov 
model yields a shock temperature of 0.461o;m> keV, an ionization timescale of 
n e t = 2.1^3 x 10 12 cm" 3 s, and half-solar abundances (0A5^qH). Assuming 
Sedov dynamics gives an average preshock H density of 1.7 ± 0.3 cm~ 3 . The dy- 
namical age estimate is 6500 ± 600 yr, while the best fit n e t value and derived 
n e gives 8200 ± 1700 yr; the weighted mean of the age estimates is 7600 ± 600 yr. 
We estimate an X-ray luminosity (0.25-4.5 keV) of (1.2 ± 0.2) x 10 37 ergss _1 (ab- 
sorbed), and (1.7 ± 0.3) x 10 37 ergss _1 (unabsorbed), in good agreement with 
the recent XMM-Newton determination. No significant excess hard emission was 
detected; the luminosity < 1.2 x 10 35 ergss -1 (2-8 keV) for any hard point source. 

Subject headings: galaxies: individual (M33) — shock waves — supernova rem- 
nants 



1. Introduction 



Multiwavelength studies of nearby galaxies are becoming increasingly effective at provid- 
ing statistically interesting samples of many classes of obiects, including supernova remnants 
(SNRs). At a distance of 817±58 kpc (Freed man et al. Il200ll ) and with a relatively face-on 
orientation {i = 55° ± 1°; IZaritsky et al. Ill989l ). the late-type Sc spiral galaxy M33 is a key 
galaxy for such studies. At the assumed distance to M33, 1" subtends 4 pc, allowing the 
morphology of some objects to be studied and many confused or crowded regions to be at 
least partially resolved, depending on the wavelength band, instrumentation, and available 
spatial resolution. 

Nearly 100 SNRs have been identified in M33, based on a combinati o n of radio and op 



tical i maging and spectroscopy (IDodorico et al. 



1985 



Viallefond et al. I Il986l ; lLong et al. 



1990 



1978 



Sabbadin 



1979; 



Smith et al. 1993; iGordon et al. 



Blair Sz Kirshner 



1998 



19991 ) . M33 has also been surveyed by each of th e imaging X-ray missions, including Ein- 



stein (lLong et al. 



Long et al. 



1996 



Misanovic et al. 



1981a 



Trinchieri et al. Ill988l ). ROSAT ( Schulman fc Bregman 



Haberl fc Pietsch IboOlh . and XMM-Newt o n JPietsch et al. I I2OO3I. 



1995 



2004 



2003). The ROSAT su rvey of lLong et al. I f !996h found 12 of the 98 op- 
tically identified SNRs in the catalog of Gordon et al. I ( 1998 . hereafter GKL98), and the 
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XMM-Newton surve y of IPietsch et al. I (12003 



20041 ) brought the total number of X-ray SNR 



identifications to 21. iGhavamian et al. I (120051 ) used archival Chandra data on M33 in com- 



parison with the optical data sets to detect X-ray counterparts to 22 of the 78 GKL98 SNRs 
within the Chandra fields of view available at that time. Using the characteristics of this 
X-ray sample, X-ray sources without optical or radio counterparts but with X-ray hardness 
ratios similar to those of confirmed SNRs were also identified as candidate SNRs. This earlier 
work was largely responsible for motivating the Chandra ACIS Survey of M33 (ChASeM33), 
a Very Large Project to exa mine the X-ray point - and extended-source populations and dif- 
fuse X-ray emission in M33 (jSasaki et al. 1 120071 ). At this writing, the survey is in progress 
and nearly complete. Here we choose to highlight observations of a single region in M33 to 
demonstrate the benefits that this survey will ultimately provide over much of the galaxy. 



M33 SNR21 (SNR #21 in the GKL98 catalog; also #121 in IPietsch et al. 1 12004 and 



#108 in iMisanovic et al. 1120061 ) is the brightest X-ray SNR in M33 and is located in the 
outskirts of the giant Hn region NGC592, which is ~9' due west of the galaxy's nucleus 
(see Fig. [1]). The relation between the SNR and the Hn region is shown in more detail in 
Figure [2, an RGB composite image constructed from continuum-subtracted Local Group 
Galaxy Survey data (LGGS3; iMassev et al. 1120061 ). using Ha (red), [Sn] (green), and [Olll] 
(blue) narrow band images; see |J3] for further details. The H II region consists of a pair 
of bright cores separated by ~ 25" which are surrounded by extensive (~ 2' x 2') faint 
filamentary structure. As indicated, the SNR stands out by vir tue of its elevated [S II] 
emission relative to Ha, a signature of shock-heated gas (see, e.g., iBlair fc Long 1 120041 ) . 



^ttp : / / www . lowell . edu/users/massey/lgsurvey . html 



-4- 




Fig. 1. — Large squares indicate the 16f8 x 16'8 ACIS-I detector footprint for the pointings 
closest to the SNR (Fiel d 4), superposed on a deep Ha image taken with the Burrell Schmidt 
telescope at Kitt Peak ( McNeil Sz Winkler 20061 ) . The small rectangle indicates the region 
covered by Figure [2] 



: 



f 



22.0 18j0 14j0 1:3G:10J) 06.0 

RA (J2000) 

Fig. 2. — False color RGB composite of continuum-subtracted images from the LGGS of the 
region surrounding M33SNR21; Ha is shown in red, [Sll] in green, and [Olll] in blue. The 
SNR appears in projection against the giant Hn region NGC592. The two bright cores of 
the H II region are seen in magenta, and some of the faint outer emission in the northeast 
can also be seen. M33SNR21 appears as the indicated cyan structure. 



Fig. 3. — Top left: Exposure-corrected (0.35-4 keV) Chandra image, binned by 0.5 ACIS 
pixels in each direction and smoothed with a 3-bin radius Gaussian. The coordinate axes 
are labeled with J2000 coordinates. The red contours are from a similarly smoothed image 
for total counts, ranging from 2-18 cts/bin in steps of 4 ct/bin; the background level is 
0.03 ct/bin. The white ellipse shows the elliptical shell model dimensions. The central 
white "+" is the X-ray center, and the white "x" to the west is the centroid position for 
the nonthermal radio source (see text). Top right: The result of 10 iterations of the Lucy- 
Richardson deconvolution algorithm. The slight shift of the bright ridge to larger radii is an 
artifact. Bottom panels: False color RGB composite of the optical images, with Ha (red), 
[Sn] (green), [Om] (blue). Bottom left: Continuum subtracted LGGS images, overlaid with 
X-ray contours. Bottom right: LGGS images prior to continuum subtraction, overlaid with 
X-ray contours. The numbers identify stars from the LGGS catalog (see text). 
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Long et al. I (ll981at ) first detected M33SNR21 in X-r ays with Einste i n, and suggested 
that it was an SNR based on its extremely soft spectrum. iGordon et al. I (119931 ) confirmed 
the SNR identification by combining radio (VLA and Westerbork), optical (imaging and 
spectral), and X-ray (ROSAT) data. The nonthermal nature of the radio emission was 
importa nt for confirming the SNR in the midst of otherwise thermal emission from the Hll 



region. iGordon et al. I (119931 ) measured the density-sensitive optical [S n] lines and found 
n e = 270 cm -3 , indicating that a dense preshock environment is likely responsible for the 
high X-ray emissivity. The high luminosity and dense environment are also an indication 
that the SNR is inde ed embedded within the Hll region and not just seen in projection 
( IGordon et al. Ill993l ). The ROSAT data allowed estimates of ~500 km s _1 for the shock 
velocity and ~4x 10 6 K for the postshock temperature. 



More recently, IPietsch et al. I (12004 ) observed M33SNR21 as part of a deep XMM- 
Newton survey of M33 and obtained an absorbed 0.2-4.5 keV X-ray flux of (1.41 ± 0.02) x 



10 13 erg s 2 cm 2 , which implies an absorbed X-ray luminosity of (1.13±0. 14) x 10 



37 —1 

erg s 



at our assumed distance to M33. The SNR has been cataloged and called by a variety 
of names, including M33X-3, 2E 0130.3+3023, 013022+30233, 022+233, GKL21, G98-21, 
GKL98 21, GDK29, RX J0133. 1+3038, and XMMU J013311. 6+303841. 

We provide below a detailed analysis of the ChASeM33 data for M33SNR21. The 
data reduction and processing steps are described in Jj2j We select the portion of the data 
with the best imaging resolution and compare with optical data for the region, finding very 
different morphologies for the optical and X-ray emitting components of this SNR (§3]). We 
then combine the high resolution X-ray data with other portions of the ChASeM33 data (in 
which the SNR was farther off-axis) to improve the statistics and perform X-ray spectral 
analyses (SHJ). We derive global SNR parameters assuming a Sedov model, and compare and 
contrast M33SNR21 with similar SNRs in the LMC - N49 and SNR 0506-68.0 - in JJ5J The 
last section, 351 summarizes our results and conclusions. 



2. X-ray Observations and Data Reduction 

The ChASeM33 survey was designed to cover the inner and most crowded regions of 
M33 with seven fields, e ach with a total exposure of 200 ks split into two observing intervals; 



see (ISasaki et al. 1 120071 ) for further details on the survey. ACIS-I is the primary detector, 
while ACIS-S2 and S3 provide additional (though far off-axis) coverage of portions of the 
galaxy adjacent to the primary positions. At this writing, both ~100 ks Field 4 pointings 
closest to the M33SNR21 position have been performed (see Fig. [TJ). The first Field 4 
pointing was split into two pieces, ObsIDs 6382 and 7226, totaling 97.2 ks after screening 
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for high background. In these data, the SNR is ~1.7' off-axis on the ACIS-I2 chip and at 
the same nominal pointing and roll (differing by less than 0.1" in both); we merged these 
datasets for subsequent analyses. The second Field 4 pointing, ObsID 6383, was not split 
and totaled 91.4 ks after screening; the SNR was ~2.1' off-axis on the ACIS-I3 chip. The 
Field 4 observations are all close enough t o on-axis to attain good spatial resolution. We 
simulated the Chandra PSF with ChaRT@ JCarter et al. Il2003h using the fitted M33SNR21 



spectrum (see Q and obtained half-power diameters of 0.93" and 1.01" for combined ObsIDs 
7226+6382 and ObsID 6383, respectively. 

M33SNR21 was also observed in five other pointings at off-axis angles ~ 8'-18' (see 
Tabled]). Although the rapid degradation of Chandra's spatial resolution at large off-axis 
angles (Chandra POG§|) makes these observations unsuitable for spatial studies of this SNR, 
the source is bright and isolated enough from other nearby sources to make these data useful 
for investigating the spectrum of the SNR as a whole. Moreover, ObsIDs 6380 (Field 3) and 
6388 (Field 7) imaged the SNR on the ACIS-S3 chip which has superior low-energy response 
and thus provides improved information about the softest X-rays. 

We reprocessed the X-ray data from the Level 1 files to remove pixel randomization and 
apply time- dependent gain changes. We screened for background flares using the ACIS-I3 
light curves except for the ACIS-S3 data, for which we used the ACIS-S3 light curves. We 
used CIAO version 3.3.0.1 and CALDB version 3.2.0. In the Field 4 ObsIDs 6282 and 7226, 
the source dithered across two ACIS 12 "bad" columns (162, 196), resulting in the loss of 
~ 13% of the counts from the SNR. These columns had been marked bad because of a 
slight excess background. However, the SNR spans only a short stretch of the columns, and 
a careful examination of the rest of the data for these columns showed that the spurious 
background would be negligible for the SNR extraction region (~ 0.3 ct for column 162, 
effectively for column 196). Accordingly, we restored those columns for the analysis to 
improve the statistics. All the Field 4 data were merged for the imaging analysis to improve 
the statistics. We appli ed a subp i xel eve nt repositioning (SER) correction to the data using 



software developed by iLi et al. I (120041 . software available from the Chandra contributed 



software page^f) to improve the spatial resolution. The SER algorithm uses the distribution 
of pulse height amplitudes (PHAs) within the 3x3 pixel event islands to improve the event 
centroiding. 

From the processed event list, we construct a 0.35-4 keV "total counts" image (ct/bin) 



^http : //cxc . harvard . edu/ chart/ 

3 The Chandra Proposers' Observatory Guide, Version 8.0 



http : / / cxc . harvard . edu/ cont-sof t/ soft-exchange . html 
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by binning the events from the merged ~ 200 ks Field 4 data set by 1/2 ACIS pixel in each 
direction (~ 1/4" x 1/4" bins). We divide this counts image by a corresponding exposure 
map constructed using the fitted SNR spectrum (see §5]) to generate an exposure corrected 
image (phot cm~ 2 s _1 ). This exposure-corrected image is then smoothed with a 3 bin radius 
Gaussian kernel. The image was overlaid by contours from a similarly smoothed version of 
the original total counts image, with contours ranging from 2 ct/bin to 18 ct/bin at intervals 
of 4 ct/bin. The corresponding background level is ~ 0.03 ct/bin. The processed X-ray 
imaging data are shown and compared with optical images of the region in Figure [3J see §3] 
for a discussion of the optical data. The upper left panel of Figure [3] shows the smoothed, 
exposure-corrected 0.35-4 keV image from the combined ~ 200 ks Field 4 data set, binned 
to 1/2 ACIS pixel in each direction, with contours from the smoothed total counts image as 
described above. The X-ray imagery reveals a strong asymmetry in the surface brightness 
distribution. The bright clump in the east has a surface brightness ~ 5 times that of the 
fainter portions of the SNR. To examine potential finer spatial structure of M33SNR21, we 
experimented with Lucy-Richardson deconvolutions of this image. The upper right panel of 
Figure [3] shows the result for 10 iterations with the CIAO tool arestore using the combined 
Field 4 (ObsIDs 6382, 7226, and 6383) total counts image, and merged ChaRT simulations 
for the PSF. The deconvolution shows an arc of bright emission in the east, and possibly 
some patchy structure within the interior. The slight shift of the bright ridge to larger radii 
is an artifact of the reconstruction. 



3. Imaging Analysis 

For optical images of the M33SN R21 field, we took data from the Local Group Galaxies 



Survey (LGGS; iMassey et al. 1 120061 ) which used the KPNO 4m telescope and Mosaic CCD 
camera to survey most of M33 (3 overlapping fields) through narrow-band Ha, [O ill] A 5007, 
and [S n] AA6716, 6731 filters, plus broadband UBVRI. M33SNR21 appears in both the north 
and central fields. From each Mosaic image we clipped out a small section centered on the 
SNR and precisely aligned these using several dozen field stars. We then selected the images 
with the best seeing in each of the 3 emission lines (all from the north field) and in the V and 
R bands (both from the central field) and matched the point-spread functions. Finally we 
scaled the broadband images appropriately, and performed continuum subtraction to remove 
the stars as effectively as possible. In this case, we subtracted V from [Olll], and R from 
both Ha and [S n] . The result is shown in Figure [2J The lower left panel of Figure [3] is an 
enlarged version of the same image with colors rescaled to emphasize the SNR, and the lower 
right panel is an RGB composite of the narrow-band images prior to continuum subtraction 
to show the locations of the bright stars. X-ray contours have been plotted on both to 



2 4 60 2 4 6 
radius (arcsec) radius (arcsec) 



Fig. 4. — Elliptical shell fit to Field 4 image data. The panels compare the data and projected 
model for the four quadrants centered on North, South, East, and West. 

facilitate comparison with the upper panels. Based on the V magnitudes and UBVRI colors 
obtained by the LGGS, the three stars seen in projection (labeled "1", "2", and "3") and the 
star furthest to the northeast ("5") are O stars in M33. Typical O star X-ray luminosities 



are < 10 33 erg s _1 (IVaiana et al. Ill98ll ) and an O star would be undetectable in the 500 ks 
of data used in our analysis. The "star" closest to the northeast rim , "4" , consists of three 
stars, two of which appear to be O stars. The third (much fainter) star has an unusual color, 
but no X-ray events were detected in the 190 ks of the Field 4 imaging data (Fig. [3]). The 
X-ray data are not significantly contaminated by stellar emission. 

To characterize the overall shape of the X-ray SNR, we fitted the 0.35-4 keV image data 
with an elliptical shell model (assuming that the symmetry axis lies in the plane of the sky), 
with density within the shell varying azimuthally as 

n(0) oc no(l + 0.5,4(1 + cos(# - 0„))). (1) 

Figure H] shows four projections of the model onto the data (including the effect of the 
simulated Chandra PSF) for quadrants centered on north, south, east, and west. The 
greater brightness in the eastern quadrant compared to the western quadrant is evident. 
The model shows the SNR to be slightly elliptical (axis ratio 1.07) with the position angle 
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Table 1. Chandra X-Ray Observations 



r leiu 


US1US 


da 
u 




ppn 


OOU.I1LS 






(arcmin) 


(ks) 




(0.35-4 keV) 


4 


7226, 6382 


1.66 


97.2 


12 


1850 




6383 


2.06 


91.4 


13 


1740 


5 


6184, 7170, 7171 


7.77 


100.1 


10 


1470 


1 


6376 


8.51 


93.1 


11 


1430 


3 


6380 


11.26 


89.7 


S3 


2155 


2 


6378 


16.21 


103.7 


S2 


1746 


7 


6388 


17.60 


88.7 


S3 


2470 



a 0fT-axis angle from the optical axis. 

b Estimated background-subtracted counts (0.35-4 keV) 




Fig. 5. — Unsmoothed total counts image, binned to 1/2 ACIS pixels. The Field 4 Bright 
Eastern ("BE") and Faint Western ("FW") wedge extraction regions are indicated. The 
extraction region for the full SNR is indicated by the ellipse labeled "SNR" . 
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Fig. 6.— Chandra ACIS-I and ACIS-S3 spectra and best fit sedov model of M33SNR21; 
the curves which are higher at low energies are ACIS-S3, and the rest are ACIS-I data. 



of the major axis ~ 9° east from north. The semimajor (semiminor) axis is 2.65"(2.48") 
with an uncertainty of ~ ±0.05". At the assumed distance of M33, this corresponds to a 
semimajor (semiminor) axis of 10.5 pc (9.8 pc) with an uncertainty of ~ ±0.2 pc, and an 
overall size uncertainty of ~ 0.7 pc when the distance uncertainty is included. The surface 
brightness peaks at position angle ~ 104° east from north. The fitted value for the fac- 
tor A in the elliptical shell model is 1.19; this implies the ambient density on the bright 
side is about a factor of 1 ± A m 2 higher than that on the faint side. The outline of 
the elliptical shell model is indicated by an ellipse in the upper panels of Figure O Com- 
pared to the fitted ellipse, the X-ray emission extends slightly further to the south than 
other directions (Fig. [3j upper left panel; see also the comparison between the northern vs. 
southern quadrants in Figure H]). The fractional thickness of the shell, AR S /R S , is poorly 
constrained because the SNR is only barely resolved; we found that a value of 0.065 provides 
an adequate description for the limb-brightening of the SNR's shell. The resulting center 
position is (a J20 oo, 5j 2 ooo) = (01 h 33 m ir!73, ±30°38'41'. / 9), indicated by the central "±" in each 
panel of Figure El This position agrees well with the XMM-Newton position determination 



flMisanovic et al. 1 120061 ) of (a J20 oo, 5 J20 oo) = (01 h 33 m ll f :76, ±30°38'42'.'l), well within the 



XMM-Newton uncertainty of 0.55". 



Using techniques similar to those in the IWeisskopf fc Hughes I (120061 ) review paper, we 



also examined the merged Field 4 hard band (3-8 keV) image for any indications of a point 
source or emission from a plerion. The extraction region for the whole SNR includes only 
14 counts in the 3-8 keV band during the entire 190 ks integration. The best fit XSPEC 
sedov model (see §1]) predicts 11 ± 5 counts in this band, so there is no evidence for a 
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significant excess of hard emission. The brightest potential point source has a total of 
3 counts (3-8 keV) within two adjacent pixels, providing a 3a count rate upper limit of 
< 4.2 x 10 -5 ct s _1 . If we assume a Crab-like powerlaw spectrum (powerlaw index a = 2.05), 
the same absorbing column as for the best fit sedov model and the response for the ACIS-I2 
observation, this count rate limit allows an unabsorbed flux upper limit for a point source 
of < 1.5 x 10~ 15 erg cm -2 s™ 1 (2-8 keV) corresponding to an unabsorbed X-ray luminosity 
upper limit of < 1.2 x 10 35 ergs _1 for the same band. 



The centroid of the nonthermal radio source ( IGordon et al. II1993I . precessed to J2000) 
is 2.8" west of the center of the X-ray emission, and indicated by a "x" toward the right in 
each panel of Figure [3] Thermal emission from the H II region may confuse the radio picture 



( IGordon et al. Ill993l ). The resolution of the radio data is ~ 7" so the ~ 3" offset of the 



radio centroid from the X-ray center is unlikely to be significant. 

The relative morphologies of the optical and X-ray emissions from M33SNR21 are in- 
teresting and somewhat perplexing. The X-ray SNR appears nearly circular, although with 
the strong brightness gradient already described, but the optical SNR appears somewhat 
more structured and extended. The bright optical emission agrees well with the outer X-ray 
contour shown in Figure [3] along the northern and northwestern X-ray limbs, apparently 
validating the relative astrometry used to align the images. However, the bright eastern 
X-ray limb is almost devoid of optical emission (although this comparison is hampered by 
the stellar emission) and faint optical emission extends farther to the southwest from the 
X-ray contour. Taken on its own, the optical morphology in the southwest is reminiscent 
of a "blow-out," which is consistent with the idea that the density is lower on the side of 
the SNR away from the bright Hn cores. The absence of optical emission on the eastern 
X-ray peak does not appear to be due to obscuration, as the iV H column derived below is 
not extreme. 
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4. Spectral Analyses 

The combined Field 4 data (the only X-ray data for which the SNR is resolved well 
enough to extract spectra from individual subregions) has ~ 3600 counts. We extracted 
spectra from two wedge-shaped regions, each containing roughly half the total counts: the 
bright eastern part ("BE", ~ 1600 counts), and the fainter western part ("FW", ~ 2000 
counts); see Figure [51 The BE region has 28% of the area. We also extracted spectra for 
the SNR as a whole using an elliptical region ("SNR"; see Figure [5]). Because ObsIDs 6382 
and 7226 fell on the ACIS-I2 chip, while ObsID 6383 fell on ACIS-I3, we extracted separate 
spectra for the merged 6382+7226 data and the 6383 data. The data were grouped to a 
minimum of 25 counts/bin, and the x 2 statistic was used for the fits. Background spectra 
were extracted from larger adjacent regions. In order to improve the fit statistics for the 
spectrum for the SNR as a whole, we also extracted spectra for a number of the far off-axis 
observations: Field 1 (ACIS-I1), and Fields 3 and 7 (ACIS-S3). The source and background 
extraction regions were enlarged appropriately to account for the expansion of the Chandra 
PSF at large off-axis angles. This provided a total of ~ 10000 counts for the SNR as a whole. 
For the ACIS-S3 data (Fields 3 and 7), the SNR was imaged near the edge of the chip and 
~ 10% of the flux spilled off the chip. 

We fitted the spectra using XSPEC version 11.3.2. In each case, we included absorp- 
tion (phabs) corre sponding to a Galactic column of Nn = 5.5 x 10 20 cm -2 toward M33 



( IStark et al. 1119951 ); this component was frozen during the fits. To account for absorption 
within M33, we included a second absorption component (vphabs) which was allowed to 
vary: the abundances of th e vphabs component were fixed to be 0.5 solar, where we assume 



Anders &: Grevesse I ( 119891 ) values for the solar a bundance set. The M33 abundance at the 
galactocentric radius of M33SNR21 is ~ 0.5 Z Q ([Henry fc Howard Ill995l ). For the thermal 
component, we examined collisional ionization equilibrium models (which provided very poor 
fits), and a number of nonequilibrium ionization (NEI) models: nei, pshock, and sedov. 
The nei and pshock models represent impulsive heating to a constant temperature, T; in 
the former model, the spectrum evaluated at a single ionization timescale, r = n e t, while in 
the latter, the spectrum is integrated over ionization timescales from r = to n e t. In the 
sedov model, the nonequilibrium ionization spectrum is integrated over a Sedov-stage SNR 
profile, and the ionization timesc ale is r = re R to where re f „,, is the postshock electron density 



and t is the age of the SNR. See iBorkowski et al. I (120011 ) for a more detailed description of 
the models. 

We performed fits for regions BE, FW, and the whole SNR using the high spatial 
resolution Field 4 data, and also for the whole SNR including data from Fields 1, 3, and 7 
in addition to Field 4. The spectra were fitted simultaneously with the model parameters 
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Table 2. Spectral Fit Results for M33SNR21. 



Parameter Units 



BE 

F4 data 



FW 
F4 data 



SNR 
F4 data 



SNR 
Fl,3,4,7 data 



const (phabs(vphabs(pshock) ) ) model 



Abundance 


(solar a ) 


n 4+ - 6 




0.3 ±0.1 


0-28±g5g 




[keV] 


n cc+o.04 

u - dd -0.07 


0.6118S 


0.58 ±0.03 


0.58 ±0.02 


r 


[10 12 cm -3 sec] 


7 +0 - 4 

u - ' -0.2 






0.6 ±0.1 


-^H,M33 b 


[10 20 cm- 2 ] 


< 7.8 


o 4+14.3 
^•^-2.0 


< 3.0 


< 2.9 


K c 


[10 -4 cm" 5 ] 


0.54 ±0.02 


fil+ - 03 

U ' U± -0.02 


1 36+ 003 

1 - ou -0.04 


1 46 +0 - 03 

± -^ u -0.02 


Xled (d.O.f.) 




1.11(51) 


1.11(62) 


1.14(111) 


1.14(259) 



const (phabs(vphabs(sedov) ) ) model 



Abundance 


(solar a ) 


5 +0 - 6 

U.d_ 2 


7+ 2 - 1 

u -'-0.3 


5 +0 - 3 

U.O_ 02 


45+ - 12 


kT e 


[keV] 


u -^°-0.05 


n A o+0.04 


0.46 ±0.02 


46 +0 - 01 

u -^°-0.02 


T 


[10 12 cm -3 sec] 


2 4 +1 - 5 

z -^-0.8 


2 4 +1 - 7 


o fi +0.9 
z -°-0.6 


9 i +0.2 
z - J --0.3 


^Vn,M33 b 


[10 20 cm" 2 ] 


< 9.6 


4 1+17.5 


< 4.0 


< 3.4 


K c 


[10 -4 cm" 5 ] 


0.51 ±0.02 


0.49 ±0.02 


1.08 ±0.03 


1.20 ±0.02 


Xled (d-O.f.) 




1.02(51) 


1.05(62) 


1.07(111) 


1.08(261) 



' lAnders fe Grevesse I (119891 ) 



b Nn,MW was fixed at 5.5 x 10 20 cm 2 . The abundances for -/Vh,m33 were fixed at 0.5 solar. 

formalization K = (10 -14 /(4tt.D 2 )) J n e n-ftdV where D is the source distance (in cm), 
n e is the electron number density (cm -3 ) and Tin is the hydrogen number density (cm -3 ). 



Note. — Error ranges are 90% confidence intervals (Ax 2 < 2.706 for one parameter). 
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(including model normalizations) linked. To account for the loss of flux off the chip for the 
ACIS-S3 data, an additional multiplicative const model was applied, allowed to be free for 
the ACIS-S3 datasets, and fixed at 1.0 for the other datasets. The resulting const fit values 
were ~ 0.9 for the ACIS-S3 data, as expected. 

All three models had comparable x 2 re d ~ 1; but the sedov models consistently produced 
better fits: 1.08 for 261 degrees of freedom (dof) for sedov vs. 1.14 (259 dof) and 1.21 (259 
dof) for pshock and nei, respectively. The nei and pshock models gave lower abundance 
estimates (0.15-0.19Z Q and O.22-O.34Z , respectively) compared to the sedov model (0.36- 
O.57Z ). The temperatures were 0.52-0.56 keV (nei), 0.56-0.6 keV (pshock), and 0.44- 
0.47 keV (sedov). The pshock and sedov fit parameters are listed in Table [2] the listed 
parameter uncertainties are 90% confidence intervals (A% 2 < 2.706 for one parameter) based 
on the statistical errors. 

In the sedov model, the temperature parameter is the postshock (electron) tempera- 
ture. In the Sedov self-similar solution the temperature increases (and the density decreases) 
radially inwa rds so that the emission -weighted X-ray temperature is ~ 1.3 times the shock 



temperature (IRappaport et al. Ill974j ). In the other two models, nei and pshock, the tem- 
perature parameter is the (constant) postshock temperature. If the temperature in the SNR 
is in fact increasing radially inwards as in the Sedov solution, the fitted temperature pa- 
rameters for the nei and pshock models would reflect an average X-ray temperature higher 
than the temperature at the shock front; this would be consistent with the generally higher 
temperature parameters in those fits compared to the sedov model fits (Table [2]). 

Based on the best fitting sedov model, we evaluated the X-ray flux and luminosity 
(absorbed and unabsorbed) for 0.35-3 keV (the range best constrained by the data) and a 
broader 0.25-4.5 keV range; the results are presented in Table |3J The quoted errors include 
the normalization uncertainty (Table [2]) and the distance uncertainty. These luminosities 
may be compared to the luminosities (sc aled to 817 kpc) e stimated with Einstein IPC 



1.7 x lO^ergs" 1 absorbed, 0.15-4.5 keV; lLong et al. Ill981ah . ROSAT PSPC (flux 1.4 x 



Table 3. X-ray fluxes and luminosities based on the sedov model fits 



Band 


F (absorbed) 


F (unabsorbed) 


Lx (absorbed) 


Lx (unabsorbed) 




(10 -13 ergs -1 cm -2 ) 


(10~ 13 ergs -1 cm -2 ) 


(10 37 ergs" 1 ) 


(10 37 ergs -1 ) 


0.35-3.0 keV 


1.46 ±0.03 


1.95 ±0.03 


1.2 ±0.2 


1.6 ±0.3 


0.245-4.5 keV 


1.50 ±0.03 


2.15 ±0.04 


1.2 ±0.2 


1.7 ±0.3 
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10 13 erg cm 2 s 



1 ; 



implying L x 



20011) and XMM-Newton EPIC (flux (1.41 ±0.02) x 10~ 13 erg cm ~ 2 s 



1.1 xlO 37 ergs -1 absorbed, 0.1-2.4 keV; lrIaberl fc Pietsch 

1 implying Lx = (1.13± 
0.16) x 10 37 ergs _1 absorbed, 0.2-4.5 keV; iPietsch et al. II2004I ). Our luminosity estimate of 
L x = (1.2± 0.2) x 10 37 erg s" 1 (absorbed, 0.25-4.5 keV) is in good agreement with the ROSAT 
and XMM-Newton estimates and in reasonable agreement with the Einstein estimate. It is 
also worth noting that th e SNR X-ray luminosity function for M33 does not reach as high as 
in other nearby galaxies (jGhavamian et al. 1 120051 ) . M33SNR21, the brightest SNR in M33, 
would be only the fourth or fifth brightest S NR if it wer e in the Large Magellanic Cloud 
(LMC), based on the luminosities reported in lLong et al. I (ll981bf ). 



Discussion 



5.1. Association with NGC 592 



As noted in (JU the SNR lies along the line of sight to the bright Hn region NGC 592. 
The high X-ray luminosity of the SNR and the relatively high density interstellar medium 
(ISM) inferred f rom the optical [S ill emission imply that the SNR is in fact embedded in 
the Hn region ([Gordon et al. II1993I ). The fact that the X-ray emission from M33SNR21 
is brightest on the side toward the bright Hn region core further supports the association 
with NGC 592. This suggests that the M33SNR21 progenitor could have been a massive 
star associated with the burst of star formation which produced NGC 592. 



Pellerin I (120061 ) recently reported on FUSE observations of NGC 592 in which the FUSE 



30" aperture covered the bright Ha cores (but excluded M33SNR21). The line profiles show 
extended blue absorption wings, indicating the presence of evolved O stars. Her fits to 
synthetic spectra based on stellar population models are consistent with a 4.0 ±0.5 Myr 
popul ation at solar metall icity; the FUV models are consistent with a ra nge of ~ 0.4-1.2 
solar. iBosch et al. I (120021 ) estimated a population age of > 4.5 Myr, but iPellerin I argues 
that this is inconsistent with the strength of the P-Cygni features. If we assume that the 



M33SNR21 progenitor belongs to the same population, the models of iMaeder I ( I1990l ) suggest 
an initial stellar mass of ~ 50 M & for half-solar to solar abundances. Our X-ray spectroscopic 
analysis of M33SNR21 with Chandra data also suggests roughly half-solar abundances. The 
lack of an X-ray signature from ejecta is likely a combination of the age of the SNR (~ 5- 
10 kyr) and the limited Chandra spatial resolution at the distance of M33. For the integrated 
spectrum for the SNR as a whole, most of the X-ray emission arises from the swept-up 
interstellar medium in M33. These high spatial resolution X-ray observations apparently do 
not provide sufficient resolution and statistics to separate out high abundance clumps if they 
are present. 
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Pellerin I (120061 ) also obtained an intrinsic extinction of E(B — V) t 





0.07 ±0.02 i 


or the 




Fitzpatrick 


1986h 


Bohlin et al. 


1978D 



implies Nn fa 1.7 x 10 21 cm -2 , while the Galactic extinction cur 
implies iVn ~ 4 x 10 20 cm -2 . The analysis of the Chandra X-ray spectra provides primarily 
upper limits Nn < (2-10) x 10 20 cm" 2 for the M33 column, consistent with a solar metallicity 
extinction curve. However, given the complex environment, the position of M33SNR21 along 
the line of sight to NGC 592 cannot be determined at present. If the SNR is toward the near 
side of the H II region, the low column for the X-ray data could also be consistent with the 
FUV data and a half-solar metallicity extinction curve. 



5.2. Sedov Analysis 



In the following, we examin e the global proper t ies o f M33SNR21 based on the as- 
sumption of a Sedov-stage SNR. iHnatyk fc Petruk I (119991 ) investigated the properties of 
nonspherical SNRs using approximate analytic models, including the cases of an SNR in 
an exponential density gradient, and an SNR in an off-center r -2 powerlaw distribution (for 
example, an explosion in the wind cavity of another star). The integrated spectral properties 
of the SNRs were found to be approximately like those of Sedov models, even though the 
detailed properties are expected to vary spatially across the SNRs. The surface brightness 
(oc n 2 ) varies more strongly with density around the remnant than the temperature (oc n), 
and the shock velocity varies only weakly (oc n 1 ^ 2 ). The properties derived from the sedov 
model fits should represent values averaged over the SNR. 

The Sedov solution depends on three parameters: the explosion energy, E , the ambient 
mass density, p , and the time since the explosion, t . The ambient density of hydrogen, n , is 
related to the mass density by n = p n /(lA3 amu), where we ass ume a fully ionized gas with 
half-solar abundances (based on the lAnders fc Grevesse I Il989l solar abundance set). The 
XSPEC sedov model fits provide estimates for the postshock temperature, T s , the ionization 
timescale, r = n e s t , and also a model normalization K = (10~ 14 /(47r.D 2 )) J dVn e nn. 
In the ionization timescale, n ejS is the postshock electron number density (where we are 
assuming electron-ion equilibration in the shock), and n e , s ~ 4.8 x no for our assumed 
abundances, and assuming strong shock jump conditions. For M33SNR21, we also have 
an estimate for the physical radius of the SNR, R s , based on the angular radius, 9 S , and 
the distance to M33. With four observa tional constrain t s and three model parameters, the 
system is overconstrained. As noted by iHughes et al. I (119981 ). the age estimates based on 



the Sedov dynamical age and the age based on ionization timescale allow for a consistency 
check. If the SNR is well modeled by a Sedov solution, these two age estimates should be 
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the consistent with each other. 

The X-ray temperature implies a shock velocity. 

1/2 



«.= (^) = 904T,'£ v kms-' (2) 

where fcT s ke v is the postshock temperature in keV and fi s is the postshock mean mass per 
free particle. The equation assumes a strong shock into a monatomic gas with three degrees 



of freedom, and complete equilibration between species (iMcKee fc Hollenbach II1980I ). For a 
fully ionized plasma with half-solar abundances, n s = (mn + A He mu e + ■ ■ ■)/(2 + 3A} ie + - ■ •) ps 
0.61, where m# and m He are the atomic weights of H and He, respectively, and A He is the 
fractional abundance of He relative to H. From the sedov model estimate for the postshock 
temperature, we have v s = 611 ± 12 km s" 1 where the error is a purely statistical error, and 
almost certainly underestimates the real uncertainty in the velocity. 

If the SNR is assumed to be expanding as R s oc t s , the shock velocity, combined with 
the current SNR radius, allows for an estimate of the SNR dynamical age, to- The model 
normalization, K, combined with v a , provides an estimate for the ambient hydrogen density, 
no, and the postshock electron density, n e , s = 4.8 n , can be estimated by using the strong 
shock jump conditions and the (half-solar) abundances. The postshock density can be com- 
bined with the ionization timescale estimate, r, from the sedov model fit to provide a second 
estimate for the SNR age. From these basic parameters, the swept up mass, Msu, and the 
explosion energy, E , can be estimated. We proceed to the evaluation of the parameters 
based on the sedov fit parameters from simultaneous fits to the data for the whole SNR (the 
right hand column in Table [2]) . 

The mean angular radius of the X-ray SNR is 2.56" ±0.13", where the uncertainty in the 
radius includes a 2% angular size uncertainty and 3.5% uncertainty (added in quadrature) 
for the azimuthal asymmetry. This angular radius and the assumed distance to M33 imply 
a radius of 10.1 ± 0.9 pc for the SNR, where the distance uncertainty is now included. For a 
Sedov-stage SNR, the deceleration parameter is 5 = 2/5, and we estimate the dynamical age 
to be (2/5)R s /v s ~ 6500 ±600 yr. Given the radial dependence of the electron and H number 
densities, the emission integral, J dV n e %, can be evaluated from the model normalization. 
For the case of the Sedov solution, 

rfVn.nH«4e(^)ngfe$) (3) 

where no is the preshock H density and £ = 0.517773 is an integration constant evaluated 
by numerical integration over the Sedov solution. The average preshock H density can be 
estimated as 

n = 1.58 K X H D 800 R~% 2 cm" 3 , (4) 
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where K_i is the XSPEC sedov model normalization in units of 10~ 4 cm -5 , -Dsoo is the 
distance to the SNR in units of 800 kpc, and R s ,io is the radius of the SNR in units of 
10 pc. From the sedov fit parameters, we obtain an average preshock H density of n = 
1.7 ± 0.3 cm~ 3 . Approximate pressure equilibrium around the rim of the SNR would imply 
that the temperature should vary as ~ tlq . Although the temperature derived from the fit 
to the BE region is somewhat higher than that obtained from the FW region fit (though not 
significantly so), even taking the full 90% confidence ranges into account does not provide 
the ~ 2:1 ratio implied by the density variation from the elliptical shell model fit (see §2]). 
The relative overpressure of the BE region compared to the FW region may be an indication 
that the encounter with the density gradient toward the east is relatively recent. This might 
also explain the lack of optical emission (see $3]) if there has not been sufficient time for 
the shocks to become radiative. In principle, the ionization timescale provides a constraint 
on the age of the encounter; however, 90% confidence intervals for r for the BE and FW 
regions are broad enough to make this a weak constraint. An additional possibility is that 
there are sufficient photoionizing photons toward the east to prevent the postshock gas from 
recombining. It is notable that the edge of the enhanced [S il] in the east coincides roughly 
with a bright rim of Ha emission associated with one of the bright cores of the H II region. 

If a Sedov-stage SNR is assumed, the ionization timescale provides for another estimate 
for the SNR age: to = T/n 6yS « 8200 ± 1700 yr. The age based on ionization timescale is 
somewhat larger than the dynamical age, but the estimates are reasonably consistent. The 
weighted mean of the estimates yields an age of 6700 ± 600 yr. 

From the preshock H density and the radius of the SNR, the swept up mass is ~ 
260 ± 80 Mq, assuming a constant preshock density (which is almost certainly not the case, 
but would provide a reasonable estimate if the density gradients are not too extreme). For 
a Sedov-stage SNR, this mass, together with the postshock temperature, implies a total 
explosion energy of E m (1.8 ± 0.3) x 10 51 erg. 



Gordon et al. I (119931 ) estimated the nonthermal pressure to be P n t/k ~ 2.4 x 10 cm K 
(where k is the Boltzmann constant) based on the radio flux and assuming a 28 pc diameter 
shell with thickness AR/R = 1/12 for the emitting region. Scaling their pressure value to 
our smaller radius and using our Sedov model par ameters to estimate t he emitting volume, 



we obtain P n t/k ~ 4.6 x 10 6 cm _3 K. As noted by lGordon et al. Ill993l . the estimated non- 
thermal pressure is much smaller than the thermal pressure estimate based on X-ray data 
(see Table H]). They note that the nonthermal pressure estimate is likely to be a lower limit: 
pressure: field/particle nonequilibrium could increase the nonthermal pressure, and volume 
for the radio emitting plasma (difficult to ascertain, because of the radio spatial resolution) 
could be much smaller than was assumed. 



-20- 



5.3. Comparisons with Other SNRs 



The asymmetric brightness distribution invites comparison with two LMC supernova 
SNRs which also show large variations in optical and X-ray surface brightness: N49, and 
0506-68.0 (also known as N23). Both SNRs show strongly asymmetric X-ray and optical sur- 
face brightnesses as well as indications of interaction with denser material. T able H] summa 



rizes the X-r ay properties and Sedov model parameters for the SNRs, based o mHughes et al. 
(jl998l . l2006h . 



Gordon et al. I (11993 ) compared M33SNR21 to the LMC SNR N49; in the latter SNR, 



both the optical and the X-ray emissio n are strongly enhanced in the southeast. A compre- 
hensive study by lVancura et al. I (119921 ) examined optical and UV data, and developed a self 
consistent model for N49. The optical emission is coming from slow (40-270 km s" 1 ) shocks 
into dense (20-940 cm -3 ) gas as the SNR interacts with a molecular cloud to the southeast, 
while the X-ray emission results from faster (730 km s -1 ) sho cks into an interclo ud medium 
(no ~ 0.9 cm -3 ). A CO cloud to t he southeast of t he SN R (IBanas et al. 1119971 ) would ac- 



count for the high densities there. IVancura et al. I (119921 ) inferred an age of 5400 yr based 
on Sedov dynamics. 

SNR 0506-68.0 also s hows strong asymmetries in its X-ray, optical, and radio morpholo- 
gies (IHughes et al. II2006I . and references therein). Such asymmetries suggest that an inter- 
action is taking place between the SNR and denser gas. However, searches for an associ- 
ated molecular cloud have not been successful; moreover, the low abs orption (as implied b y 



the low column dens ity) argues against the presence of su ch a cloud (jHughes et al. 1119981 ). 



Hughes et al. I (120061 ) noted that the open cluster HS114 (IHodge fc Sexton 1119661 ) lies near 
the high density side of the SNR (only 2 pc in projection, assuming a distance of 50 kpc to 
the LMC). They found at least a factor of 10 variation in the ambient density around the 
rim of the SNR, and also detect regions with enhanced abundances in parts of the SNR. By 
assuming that the side of the SNR which is expan ding into a regi o n of low density can be 
modeled approximately as a free Sedov expansion, IHughes et al. I (120061 ) estimated an age 
of ~ 4600 yr. 



Table 4. Comparison of M33SNR21with LMC remnants N49 and 0506-68.0. 



SNR 


R s 






n 


v s h 


P s /k c 


to 


E 






(PC) 


(10 36 ergs" 1 ) 


(10 6 K) 


(cm~ 3 ) 


( km s" 1 ) 


(10 7 cm~ 3 K) 


(10 3 yr) 


(10 51 erg) 


(M & ) 


M33SNR21 


10 


12 


5.3(0.2) 


1.7 


610 


8.3 


6.7 


1.8 


260 


N49 d 


8.2 


6.3 


6.7(0.1) 


2.6 


690 


16 


4.4 


1.5 


210 


0506-68.0 d 


6.7 


2.5 


6.2(0.1) 


1.6 


660 


9.1 


3.8 


0.46 


70 


0506-68.0 6 


12 






0.25 






4.6 







a 0.5-5.0 keV 
b From Eq. [2] 

c P s /k = 4 x (2.3no)T s , where the 4 accounts for the strong shock density jump, and 2.3 is the number of free particles per 
H in a fully ionized plasma 

^Hughes et al~l (|l998h 



Hughes et al. I (120061 ) 
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M33SNR21 shows similar asymmetries in the X-ray emission, with the rim toward the 
H II region core being much brighter than the rest of the SNR. Despite the X-ray similari- 
ties with the other SNRs, the optical emission shows rather different properties. N49 and 
SNR 0506-68.0 both show optical brightening on the same side as the X-ray brightening. In 
contrast, M33SNR21 does not show strong optical brightening near the X-ray bright eastern 
rim (though the situation is confused by stellar contamination), but instead shows an optical 
extension toward the southwest well beyond the faint X-ray rim. The strong X-ray emission 
ahead of the optical emission in the east is puzzling. As noted in §5.21 this may indicate 
a relatively recent encounter with denser gas for which photoionization associated with the 
bright H II region core maintains the ionization at a high enough state that any [S il] emission 
is weak. The extension of the optical emission well beyond the X-ray rim in the southwest 
is suggestive of some sort of breakout into a lower density region which has subsequently 
driven slower (hence X-ray faint) radiative shocks into denser material. In that case it is 
puzzling that no X-ray emission is seen within the breakout. Possibly the X-ray emitting 
plasma driving the optical shocks is tenuous enough that the X-ray emission is too faint to 
be seen even in this 190 ks exposure. In any case, the peculiarities of the X-ray and optical 
emission for M33SNR21 suggest an explosion in a complex environment. 



6. Conclusions 

We have presented the results of an analysis of the C hASeM33 project data for M33SNR 21 



an SNR embedded in the giant Hn region NGC592 (IGordon et al. Ill993l ). These data 
provide the first well- resolved X-ray imagery of an SNR in M33. The remnant is slightly 
elliptical, and roughly 5" in diameter. We fitted a three dimensional elliptical shell model 
to the X-ray distribution, from which we measured the X-ray center of the SNR to be at 
RA=01 h 33 m ll!73, Dec=+30°38'41 / /9 (J2000), and the dimensions of the SNR shell in pro- 
jection to be ~ 21.0 pc x 19.6 pc at the assumed distance to M33. The Chandra data show 
the emission to be asymmetric, with the eastern rim (closest to the bright H II region cores) 
roughly 5 times brighter than the rest of the SNR. This asymmetry suggests that the SNR 
is interacting with the H II region and further supports the argument that t he SNR is actu 



ally e mbedded in, rather than merely in projection against, the Hn region (IGordon et al. 



19931 ). The association with NGC592 suggests that the progenitor could have been a high 



mass star, and that the SNR was the result of a core collapse supernova. 

We augmented the high resolution X-ray data with other ChASeM33 data taken further 
off-axis in order to improve the statistics for X-ray spectral analyses. By fitting with the 
XSPEC sedov model, we obtained an estimate for the shock temperature kT = 0.461q 02 keV, 
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an average preshock ISM (hydrogen) density of uq w 1.7 ± 0.3 cm -3 , and an ionization 
timescale r = n es t « 2.11q 3 x 10 12 cm _3 s. Although the imaging data show the SNR 
to be somewhat as ymmetric, the global ave rage properties can be described approximately 



by a Sedov model (IHnatyk &: Petruk II1999I ). The detailed properties would be expected to 
vary azimuthally around the SNR because of the variation in the preshock density. The 
surface brightness will vary most strongly with density (oc n 2 ), the postshock temperature 
varies oc n, and the shock velocity is relatively insensitive to the ambient density (oc n 1 / 2 ). 
However, fits to the high spatial resolution data for sectors containing the brightest emission 
(region BE) and the rest of the SNR (region FW) did not show any significant differences in 
spectral properties. 

From the sedov model estimates, we can also derive average values for other properties. 
The expansion velocity of the SNR is v s ~ 611 ± 12 km s _1 based on the observed size 
of the remnant and the postshock temperature. By assuming a Sedov expansion rate, we 
obtain a dynamical age estimate of 6500 ± 600 yr. The average postshock electron density, 
n e , s = 4.8 no, where uq = 1.7 ± 0.3, combined with the ionization timescale, r, yields an 
ionization age estimate of 8200 ± 1700 yr. The estimates are reasonably consistent with each 
other, and the weighted mean of the ages is 6700 ± 600 yr. The large amount of swept up 
mass, 260 ± 80 M , also supports the argument that M33SNR21 is at least a middle aged 
SNR. 

Our fitted abundance value, 0.45loog solar, is consistent with the expected half-solar 
abundances appropriate to the galactocentric radius of M33SNR21. Although clumps of 
enhanced eject a hay e been found in other SNRs of a similar age (e.g., SNR 0506-68.0 (N23); 



Hughes et al. 1120061 ). we do not find evidence here for enhanced abundances due to ejecta. 
Fits to the SNR as a whole will dilute any ejecta signature in the much larger mass of swept 
up ISM, while much smaller extraction regions, BE and FW, yielded poor counting statistics 
for the available 190 ks of high resolution X-ray imaging data. 

The total X-ray luminosity (0.25-4.5 keV) of the source is estimated to be (1.2 ±0.2) x 
10 37 ergs _1 (absorbed), or (1.7± 0.3) x 10 37 ergs _1 (unabsorbed), in good agreement with the 
X MM-Newton re s ult of (1.13 ± 0.14) x 10 37 erg s -1 (absorbed), based on the observations 



of iPietsch et al. I (120041 ) . We searched for any emission from a point source or plerion. We 
found no evidence for a significant excess of hard emission, and obtained an upper limit of 
Lx < 1-2 x 10 35 ergs _1 (2-8 keV) for the unabsorbed luminosity of any hard pointlike source 
if present. 

We compared the X-ray images with available optical data and found some puzzling 
morphological differences. In contrast to the slightly elliptical X-ray remnant, the optical 
remnant is more oblong. The bright optical emission agrees with the X-ray contours in the 
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north and northwest, but falls off more rapidly in the east, with the region of brightest 
X-ray emission almost devoid of optical emission (although this comparison is hampered 
by the presence of several bright stars). To the southwest, the optical emission extends 
considerably beyond the X-ray emission. The lack of optical emission to accompany the 
bright X-ray emission in the east could be the result of a relatively recent interaction with 
denser material, in which the shocked gas has not had time to become radiative; the apparent 
overpressure in the BE region compared to the FW region would support that interpretation. 
The edge of the [S il] emission in the east coincides with the rim of an H II loop, and the 
bright X-ray emission may be the result of an interaction with a bubble driven by the stars 
associated with the southern Hll region core. In that case, the lack of [Sn] emission may 
be the result of the photoionizing flux from the stars in the Hll region core preventing 
the postshock gas from recombining. On the other hand, the lack of X-ray emission to 
accompany the optical emission in the southwest is also puzzling. It may be that even with 
190 ks of data, the X-ray emission in the southwest is too faint to be picked up. The complex 
multiwavelength morphology suggests that the shell is expanding into a non-uniform density 
medium, as may be expected for an SNR evolving in an active star forming region. 

This work has made use of SAOImage DS£@, developed by the Smithsonian Astrophys- 
ical Observatory (Joye & Mandel 2003), the XSPEC0 spectral fitting package (Arnaud 1996V, 
the FUNTOOLSZl utilities package, the HEASARC FTOOL£§ package, and the CIAcfl 
(Chandra Interactive Analysis of Observations) package. Support for this work was pro- 
vided by the National Aeronautics and Space Administration through Chandra Award Num- 
ber G06-7073A issued by the Chandra X-ray Observatory Center, which is operated by the 
Smithsonian Astrophysical Observatory for and on behalf of the National Aeronautics Space 
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